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摘要 : 芯 位 重建 是 地 面 宇宙 线 探测 实验 开展 物理 分 析 的 前 提 和 基础 。 本 文 基 于 广 延 大 气 往 射 


的 对 称 性 ， 结 合 探测 器 的 特征 ,提出 一 种 新 的 芯 位 重建 算法 椭圆 拟 合法 ， 用 于 模拟 数据 


的 芯 位 重建 并 与 重心 法 对 比 。 对 于 能 量 大 于 1 TeV 的 芯 内 宇宙 线 质 子 事例 ， 芯 位 分 辩 率 在 小 


于 5m, 明 显 优 于 重心 法 ; 对 于 临近 阵列 边界 20 m 范围 内 能 量 大 于 1 TeV HS MED HR 
小 于 10 m， 相 比 于 传统 方法 有 显著 优势 。 由 于 该 方法 能 够 在 保证 精度 的 前 提 下 重建 边沿 事 


例 和 部 分 芯 外 事例 ， 极 大 地 增 大 了 事例 利用 率 ， 对 于 1-10 PeV 质子 事例 的 探测 器 阵列 面 
利用 率 比 重心 法 增 大 1 倍 左右 。 
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宇宙 线 的 能 量 跨越 10 多 个 量 级 、 流 强 跨 越 30 多 个 量 级 ， 近似 符合 UN/IdEE”, 其 
中 谱 指数 y 的 取 值 约 为 2.7 叫 。 从 MeV 到 ~400 TeV 能 区 字 宙 线 流 强 较 大 ， 可 以 用 卫星 直接 
探测 ， 更 高 能 区 宇宙 线 由 于 流 强 太 低 主要 利用 地 面 实验 间接 探测 由 。 水 切 伦 科 夫 探测 器 是 
地 面 宇宙 线 探 测 技术 的 一 种 ， 它 具有 全 天 候 、 宽 视 场 等 优势 ， 主 要 是 通过 测量 原初 宇宙 线 进 
入 大 气 后 产生 的 广 延 大 气 簇 射 重建 出 原初 宇宙 线 信息 W”" 。 位 于 墨西哥 的 HAWC 和 位 于 四 
川 稻城 的 LHAASO-WCDA 等 多 个 宇宙 线 探测 实验 都 采用 了 水 切 伦 科 夫 探测 器 技术 。 

T EKAP (Extensive Air Shower, EAS) E jai fe ES AEA KAU SKA 


份 发 生 相 互 作 用 , 通过 强 子 级 联 和 电磁 级 联 产生 大 量 次 级 粒子 的 过 程 , 通常 用 纵向 发 展 和 横 
iğ 向 分 布 来 描述 其 物理 特征 。 常用 于 描述 纵向 发 展 的 模型 有 Gaisser-Hillas function °°"), 
Greisen function’**"#11 Gaussian-in-Age5?， 三 者 都 与 实验 测量 值 吻合 得 很 好 。 EAS 的 横 
= 向 分 布 则 通常 用 横向 分 布 函数 ( Lateral distribution function, LDF ) 和 
= Nishimura-Kamata-Greisen (NKG) 函数 来 描述 ps1630。 LDF 是 简单 导出 的 、 纯 数学 拟 合 的 、 
昔 述 平均 分 布 的 经 验 公 式 ， 它 的 参数 不 直接 关联 广 延 大 气 复 射 的 物理 特性 “ ”0。 NKG É 
Al 数 比 LDF 更 接近 物理 本 质 ,一 方面 NKG 函数 能 够 反映 EAS 的 横向 分 布 和 纵向 发 展 等 物理 特 


N 
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征 ， 另 一 方面 在 拟 合 过 程 中 还 充分 考虑 了 EAS 的 统计 涨 落 5 OS, 
地 面 宇宙 线 探测 实验 通过 探测 包括 电子 (e+、e-) . BF Cu. pt) 和 光子 在 内 的 EAS 
次 级 粒子 密度 分 布 和 前 锋面 到 达 时 间 反 演 事 例 的 芯 位 、 方 向 和 能 量 等 信息 ， 这 个 过 程 称 
为 宇宙 线 事例 重建 或 EAS 重建 ， 心 位 重建 是 EAS 重建 的 重要 内 容 。 常 用 的 芯 位 重建 方法 有 
© 重心 法 、 树 形 分 析 法 、 最 大 似 然 法 、 二 维 高 斯 拟 合法 等 ， 这 些 方法 各 自 都 有 独特 的 优势 和 不 
足 W% 1。 重心 法 和 树 形 分 析 法 的 优势 在 于 计算 量 小 ， 节 省 计算 资源 ;其 不 足 在 于 重建 时 
破坏 了 事例 的 对 称 性 特征 。 芯 位 重建 的 精度 依赖 重建 出 的 芯 位 到 探测 器 中 心 的 距离 。 最 大 
似 然 法 重建 芯 位 精度 高 ， 但 是 精度 也 依赖 探测 器 单元 触发 数量 Nai) ， 当 Ni 大 于 1000 时 
其 蕊 位 分 辨 率 小 于 2 米 , 精度 相 比 前 两 者 显著 提高 ; 但 是 一 方面 该 方法 计算 量 较 大 需要 占用 
大 量 计算 资源 ， 另 一 方面 其 精度 依赖 于 能 量 、 方 向 和 Noa 等 诸多 因素 。 二 维 高 斯 拟 合法 是 
一 种 简单 的 数学 拟 合 方法 ,克服 了 重心 法 和 树 状 分 析 法 的 部 分 缺点 ， 对 于 世 位 位 于 探测 器 阵 
列 内 部 的 事例 ( 蕊 内 事例 ) 重建 效果 较 好 ， 而 对 于 蕊 位 位 于 探测 器 阵列 外 的 事例 不 能 准确 重 
RECTION, 除 最 大 似 然 法 外 ， 这 些 重建 方法 都 需要 通过 筛选 条 件 丢 弃 大 量 数 据 ， 只 保留 芯 
位 在 阵列 中 心 附 近 小 范围 的 事例 ， 以 牺牲 探测 器 有 效 面 积 的 方式 来 保证 重建 结果 的 可 靠 性 ， 
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造成 大 量 观测 数据 的 浪费 5 所。 此 外 水 切 伦 科 夫 探测 器 数据 的 低 信 噪 比 也 是 重建 的 挑战 之 
一 , 如 LHAASO-WCDA 的 观测 数据 在 去 噪 之 前 其 信 噪 比 约 为 70/218, 通过 遍历 法 去 噪 之 后 
噪声 比 约 为 22/70， 也 就 是 说 对 于 触发 70 个 探测 器 单元 的 事例 ， 其 中 约 有 22 个 通道 的 信 
号 是 噪声 ,这些 噪声 会 给 精确 重建 带 来 很 大 的 困难 ， 而 上 述 方法 几乎 都 不 能 很 好 地 处 理 噪声 
带 来 的 影响 "9 。 

本 文 基于 广 延 大 气 簇 射 的 对 称 性 和 水 切 伦 科 夫 探测 器 阵列 的 物理 特征 , 提出 了 一 种 新 的 
蕊 位 重建 方法 一 一 椭圆 拟 合法 。 该 方法 充分 利用 EAS 的 对 称 性 ， 消 除 噪声 对 芯 位 重建 的 影 
响 。 在 满足 芯 位 重建 精度 的 前 提 下 ,能 够 重建 出 芯 位 落 在 探测 器 边沿 的 事例 和 蕊 外 事例 ， 扩 
大 探测 器 事例 利用 率 。 
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1 .水 切 伦 科 夫 探 测 器 技术 

水 切 伦 科 夫 光 探测 器 是 利用 纯 水 作为 介质 ， 通过 光电 倍增 管 收 集 带 电 粒 子 在 水 中 运动 
产生 的 切 伦 科 夫 光子 来 测量 广 延 大 气 簇 射 的 地 面 粒子 探测 技术 。 现 有 各 实验 采用 的 水 切 伦 
科 夫 探测 器 阵列 的 单元 探测 器 构造 和 探测 原理 基本 相同 ， 探测 器 单元 的 各 种 性 能 参数 和 布 
局 上 存在 一 定 的 差别 。LHAASO-WCDA 是 最 具 代 表 性 的 水 切 伦 科 夫 探测 器 阵列 之 一 。 当 
FERRATE WCDA 视 场 内 产生 广 延 大 气 簇 射 ， 次 级 粒子 会 进入 阵列 内 的 纯 水 中 做 超 光速 
《介质 中 的 光速 ) 运动 激发 切 伦 科 夫 光 ， 每 个 基本 单元 的 光电 倍增 管 PMT) 收集 到 的 切 
伦 科 夫 光 电子 数 (NPE) 与 该 位 置 的 EAS 次 级 粒子 密度 成 正比 已 ;探测 器 单元 触发 时 间 可 
精确 到 纳 秒 量 级 ,通过 这 些 数据 可 以 重建 营 位 、 方 向 等 信息 4。 其 他 实验 的 水 切 伦 科 夫 
探测 器 阵列 如 HAWC 的 结构 和 布局 与 LHAASO-WCDA 有 所 不 同 ,但 是 其 探测 基本 一 致 。 
LHAASO-WCDA 由 3 个 水 池 组 成 , 具有 3120 个 单元 探测 器 以 及 后 级 的 电子 学 、 定 时 、 数 
据 采 集 、 触 发 判 选 、 数 据 处 理 、 标 定 等 功能 系统 。 单 元 探测 器 (图 1) 为 5m x5 m 的 水 域 , 深 
FER 44 m， 两 个 单元 之 间 用 隔 光 帘 隔 开 ,主要 目的 是 避免 来 自 同一 次 级 粒子 ,尤其 是 绢 子 
号 的 串扰 。 WCDA 共 采 用 3120 支 大 尺寸 光电 倍增 管 PMT) 分 别 布设 于 每 个 单元 的 中 央 ， 置 
FKE, 向 上 观测 。 另外 还 在 水 池 中 放置 3120 支 小 尺寸 PMT， 小 尺寸 PMT 放置 在 大 尺寸 
PMT 旁边 , 用 于 扩大 艇 射 粒子 数 测量 动态 范围 ,从 而 实现 高 能 宇宙 线 的 高 精度 测量 1。 
本 文 将 以 150 mx150 m 的 LHAASO-WCDAI1 水 且 伦 科 夫 探测 器 阵列 为 样本 ， 阅 述 椭圆 拟 合 
法 重建 芯 位 的 物理 思路 和 算法 设计 。 
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2 . 椭圆 拟 合法 


5 m 5m 


图 1: LHAASO-WCDA 单元 探测 器 结构 示意 图 (E) 和 探测 示意 图 (A) 


Figure 1: Schematic diagram ofLHAASO-WCDA unit detector structure (left) and detection diagram (right) 
如 图 1( 右 〉 Bras, SEAS 以 一 定 的 天 项 角 入 射 到 地 面 ， 到 达 探 测 器 阵列 平面 EAS 次 
级 粒子 集中 在 很 薄 的 前 锋面 内 ， 次 级 粒子 横向 分 布 用 可 以 NKG 函数 描述 


N 工 (4.5 一 Pe 
s ( s) Cy His 4.5 
ra 2A (S) (4.5258) ry i, 


(1) p(7)= 


C) 式 中 入。 是 EAS 的 总 粒子 数 ， 主 要 依赖 于 原初 宇宙 线 粒 子 的 能 量 ，rm ER EAB, 
与 原初 宇宙 线 的 成 分 有 关 ; s 是 描述 EAS 发 展 阶段 的 参数 ， 被 称 为 EAS 的 年 龄 ， 当 s=7 时 


I(4.5— 
EAS 发 展 到 极 大 ， Fore oy 是 归 一 化 系数 。 (1) 式 表明 的 是 年 龄 为 s 的 EAS 前 


锋面 次 级 粒子 随 距 离 前 锋面 距离 + 增 大 而 减 小 , 年 龄 越 大 减 小 越 慢 次 级 粒子 密度 分 布 函数 越 
平坦 。 当 一 个 天 顶 角 为 9 的 EAS 前 锋面 接近 探测 器 阵列 时 ， 近 似 认 为 在 前 锋面 到 达 探 测 器 
的 几 十 纳 秒 内 次 级 粒子 密度 变化 忽略 不 计 ， 那 么 到 达 任 意 位 置 y) 处 探测 器 单元 的 次 级 
粒子 密度 为 : 


o(x, y) = p(r)cos o (2) 


其 中 


(x-x 7 A-sin? pcos’ 0)+(y— yp)? (l -sin? psin® 6) — 2(x—xy)(y — yo) sin’ psin 8 cos 8 
Yn tO- y) 


(3) 


(xoyo) 是 地 平面 上 的 芯 位 ，〈0, 9) 是 EAS 的 方向 ， 其 中 0 是 方位 角 , 9 是 天 顶 角 ， Coy) 


是 平面 上 任意 位 置 探测 器 单元 的 坐标 。 从 表达 式 可 以 看 出 o(X, y) EAS 的 芯 位 和 方向 有 


KR 其 数学 结构 十 分 复杂 不 方便 分 析 其 特征 ,为 了 更 加 直观 地 理解 前 锋面 次 级 粒子 分 布 情况 ， 
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进行 了 数值 计算 并 得 到 不 同 天 顶 角 和 方位 角 事例 投射 到 地 面 的 次 级 粒子 分 布 情况 如 图 2 所 
示 。 从 图 2 可 以 看 出 EAS 次 级 粒子 密度 呈现 出 明显 地 对 称 性 ， 当 天 项 角 为 0° 时 次 级 粒子 密 
度 等 高 线 为 同心 圆 , 当天 顶 角 不 为 零 时 其 等 高 线 为 同心 椭圆 , 天 顶 角 越 大 椭圆 的 离心 率 越 大 。 
因此 可 以 用 椭圆 系 方程 来 描述 次 级 粒子 密度 等 高 线 : 


[(x, —x0)cosO+(y, — yo0)sin 0J a [(x -xsin + (y, —y0)cosO _ 
2 2 = 
a, b, 


其 中 (xo,yo) 是 这 一 系列 椭圆 的 共同 中 心 ，aj,b, 分 别 是 第 n AGL TEAS ee, @ 是 EAS 
的 方位 角 ， 椭 圆 长 短 轴 之 比 awbn=fty,s) 是 EAS 天 顶 角 @ 和 年 龄 s 的 函数 。 
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图 2: 利用 NKG 函数 数值 计算 得 到 的 不 同 天 项 角 EAS 带电 次 级 粒子 密度 ”其 方位 角 都 是 45”, 其 中 a, b, c 的 天 项 角 分 别 为 0”， 
30”，45”， 色 度 图 为 次 级 粒子 密度 的 常用 对 数 ，d，e，f 分 别 为 a，b，c 的 等 高 线 表示 。 

Figure 2: The azimuth angles of the charged secondary particle density of EAS at different zenith angles calculated by using the NKG 
function are all 45°, where the zenith angles of a, b, and c are 0°, 30°, and 45°, respectively. The degree diagram is the common logarithm 


of the secondary particle density, and d, e, and f are represented by the contour lines of a, b, and c, respectively. 
WCDA 各 基本 单元 PMT 探测 到 的 NPE 与 该 位 置 的 次 级 粒子 密度 之 间 是 线性 关系 ， 
此 其 NPE 高 线 可 近似 为 椭 贺 。 用 椭圆 方程 拟 合 NPE 的 等 高 线 可 以 得 到 一 系列 的 同心 椭圆 ， 
这 些 椭 圆 的 共同 中 心 就 是 EAS 的 蕊 位 ,椭圆 的 长 短 轴 之 比 反 映 EAS 天 项 角 信 息 ， 椭 圆 的 长 
轴 反 映 了 EAS 的 方位 角 信 息 。 拟 合 NPE 的 等 高 线 系列 来 重建 芯 位 ， 一 方面 对 于 攻 位 位 于 
探测 器 边沿 使 得 EAS 对 称 性 严重 破坏 的 事例 ， 可 以 通过 局 部 信息 来 “ 补 齐 ” 顶 圆 等 高 线 〈 如 
3 所 示 ) ， 利 用 对 称 性 来 重建 出 准确 的 芯 位 ， 男 一 方面 利用 所 有 椭圆 应 该 “同心 ”的 特点 ， 
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舍弃 掉 信 噪 比 较 低 的 数据 组 ， 从 而 实现 抑制 噪声 的 效果 ， 有 具体 算法 流程 如 图 4 所 示 : 
第 一 步 : 读 取 事 例 获 取 对 心 位 重建 有 用 的 单元 探测 器 坐标 yd ， 
第 二 步 : 把 数据 按照 Npr 降序 排列 ， 并 把 数据 分 成 n 组 ， 


第 二 步 
第 三 步 : 拟 合 得 到 mn 个 椭圆 方程 得 到 nm 个 oyo) > 


»|¥o— Yor» 


: RPI Xo Fy WEZ ns Yos MASE — xo, 


第 四 步 
第 五 步 : 判断 是 否 满足 | 一 ,| < de, & Vy — Vo;|< deur 


如 果 满 足 进入 下 一 步 , R 


则 去 掉 偏 差 值 最 大 的 点 ， 回 到 第 四 步 ， 直 到 满足 条 件 或 者 n<3; 
第 六 步 : 计算 必 位 30 米内 的 贺 形 区 域 与 30 到 40 米 区 域 环形 区 域 的 Npg 总 和 之 比 NBN, 


并 且 判 断 NBN 是 否 大 于 NBNjn， 若 大 于 则 输出 上 一 步 得 到 的 (N, 3) 为 重建 蕊 位 ， 否 则 判 


断 重 建 失败 ， 该 事例 不 纳入 统计 。 


= 


图 3: 椭圆 拟 合法 重建 芯 位 示意 


Figure 3: Schematic diagram of core position reconstruction by ellipse fitting method 

该 方法 中 涉及 到 的 几 个 参数 的 选择 和 使 用 : 

1. NPE 最 大 的 点 。 对 于 忌 内 事例 ，NPE 最 大 的 点 不 能 用 于 拟 合 椭 圆 。 这 个 点 很 可 能 位 
于 椭圆 中 心 附近 。 

2. 椭圆 系 中 拟 合 椭圆 所 用 的 数据 点 的 个 数 。 离 蕊 位 越 远 的 椭圆 等 高 线 上 的 探测 器 单元 越 
多 ,因此 用 于 拟 合 一 个 椭圆 的 点 就 应 该 越 多 ,并 且 理 论 上 与 离 蕊 位 距离 成 正比 ， 因 此 我 们 采 
用 从 内 到 外 线性 增加 的 分 组 方式 来 对 数据 进行 分 组 上 文 的 推导 表明 拟 合 平面 上 任意 椭圆 最 
少 需要 5 个 点 ， 故 取 每 组 数据 为 5n (n=1,2,3-) ， 即 第 一 组 数据 5 个 点 ， 第 二 组 数据 10 个 
点 ， 第 三 组 15 个 点 以 此 类 推 。 由 此 可 得 触发 探测 器 单元 数 Ni 与 拟 合 椭圆 数 n 有 如 下 关系 : 


5n(n+1) éN,, 
2 (5) 
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Ean 
把 数据 按照 NPE 大 小 分 成 组 


个 事例 ， 把 x;,y;,NPE; 存 入 数组 


拟 合 得 到 n 个 椭圆 中 心 (Xoiyoi) » 1=0,1,2... 


TE SEF EYER). Vos PME AEX osp Vosn WIA EE, — Xo, |5|%o 一 yoi| 


3. decns MAW AEWA MT EU NPE 等 高 
由 于 所 有 的 椭圆 等 高 线 应 该 是 “同心 ”的 ， 

但 是 为 了 减 小 统计 误差 提高 蕊 位 重建 的 精度 , 需要 充分 利用 数据 拟 合 尽 可 能 
Al, 拟 合 得 到 的 椭圆 系 列 中 心 坐 标 会 有 一 定 
(1) 椭圆 系 中 心 坐 标 分 布 在 很 小 的 
坐标 总 体 密 集 分 布 、 
E 心 作为 蕊 位 《如 图 5 Ata) ; 
1 处 理 。dec 是 用 来 量化 判断 离散 度 的 标 


宙 线 事例 的 芯 位 。 
可 以 得 到 站 位 ， 


多 的 椭圆 。 由 于 物理 涨 落 和 探测 器 单 元 差异 等 原 
通常 会 出 现 以 下 三 种 情况 : 

些 点 的 重心 作为 最 终 的 重建 芯 位 坐标 ; 
需要 去 除 掉 偏离 较 远 的 那些 坐标 值 ， 
圆 系 中 心 坐标 分 布 极为 离散 ， 


的 偏差 。 


[是 ] 


Meeks 7) 


Es 
Xosp < Alou & Yosp < deny 


bta 


LE 


K Piili 


图 4: 椭圆 拟 合法 算法 流程 图 
Figure 4: Algorithm flow chart of ellipse fitting method 


准 ， 通 常 选择 deon 的 值 越 小 得 到 的 芯 位 精度 越 高 ， 


效 数据 导致 重建 效率 很 低 ， 通 
术 第 (2) 种 情况 中 ， 
HF NPE 比较 小 的 探测 器 单元 信 


重建 效率 。 在 上 述 


与 局 位 偏差 较 大 ， 这 是 


E 


起 


通常 中 低 和 


数据 组 拟 合 的 椭圆 中 心 4 
经 饱和 。 


4. xoi 和 yoi 的 标准 差 。 用 xoi 和 yo 的 标准 差 作 为 评判 参数 来 判断 这 一 


(2) 椭圆 系 中 心 
以 密集 点 群 的 
这 种 事例 作为 无 效 事 斧 


上 只 需要 所 


大 的 xoiyoi 
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线 得 到 椭圆 中 心 (xoyo) 作为 宇 


上 只 需要 


晶 是 选择 太 小 的 decu 


号 信 噪 比 低 ; 而 


坐标 与 世人 位 偏差 较 大 , 可 能 是 


拟 合 一 条 等 高 线 就 


区 域内 ， 可 以 取 这 
个 别 偏离 较 远 ， 
(3) Hi 


可 能 会 丢掉 部 分 有 


过 多 次 尝试 比 对 ， 选 择 了 Sm 这 一 数值 ， 能 够 兼顾 芯 位 精度 与 
事例 的 小 NPE 数据 组 拟 合 的 椭圆 中 心 华 标 
高 能 事例 的 大 NPE 
因为 高 能 事例 芯 位 附近 的 光电 倍增 管 已 


系列 椭圆 中 心 的 离 


散 程度 ， 若 其 标准 差 很 大 (如 上 的 (2〉 G) 两 种 情况 ) ， 说 明 一 些 数 据 组 几乎 不 包含 芯 位 


信息 ， 这 些 NPE 主要 是 由 噪声 造成 的 ， 


因此 把 这 些 数据 当做 噪声 去 除 掉 ， 而 椭圆 中 心 坐 标 


很 集中 则 表明 与 之 对 应 的 数据 中 噪声 比 很 低 , 保留 其 对 世 位 坐标 的 贡献 。 图 5 是 一 个 事例 分 
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组 重建 出 的 各 个 椭圆 中 心 坐标 ， 最 终 只 保留 红 框 内 的 密集 点 求 平均 作为 最 终 的 蕊 位 重建 结 
果 。 


x/m 


x/m 


图 5: 椭圆 拟 合法 重建 出 的 各 个 椭圆 中 心 和 真实 芯 位 ， 其 中 蓝 色 “+” 是 各 组 数据 重建 得 到 的 椭 贺 中心， 红色“*” 是 真实 芯 位 。 


Figure 5: The center of each ellipse and the real core position reconstructed by the ellipse fitting method, where the blue "+" is the center 


of the ellipse reconstructed from each group of data, and the red "*" is the real core position. 
5..NBN 与 NBNwn。 得 到 满足 离散 度 要 求 的 椭圆 中 心 坐 标 平均 值 (X60,) 后 ， 还 需要 利用 
附近 10 米 圆 和 10 到 20 米 的 环 内 NPE 总 和 之 比 NBN=NioV2o 判 断 该 点 周围 的 NPE 分 布 


\ 


符合 心 位 附近 的 NPE DER, NBN>NBN pin 2e A ES IRE DF 10m, 可 以 把 计算 结 


作为 重建 芯 位 的 最 终结 果 ， 其 中 NBNoam 的 取 值 选择 可 以 通过 数值 估计 和 数据 统计 得 到 。 在 
需要 重建 阵列 边界 附近 的 芯 外 事例 的 时 候 ， 考 虑 到 Nio 可 能 为 0， 可 以 用 NBN=NzoVio 或 
NBN=NioVz 也 能 达到 同样 的 判断 效果 ， 只 是 涨 落 变 大 且 精 度 变 差 ， 只 能 分 别 判断 重建 芯 位 
误差 是 否 低 于 20 m,30 m. 


3 . 椭圆 拟 合法 重建 结果 


蒙特 卡 罗 模 拟 是 研究 宇宙 线 物理 的 重要 手段 ， 本 文 借助 Corsika 模拟 的 广 延 大 气 簇 射 和 


Geant4 模拟 的 探测 器 响应 数据 来 验证 椭圆 拟 合法 重建 蕊 位 的 实际 效果 。 选 取 能 量 为 10TeV 
到 10PeV 的 质子 事例 ， 其 天 顶 角 为 0-40° 服 从 cos’ 分 布 ， 方 位 角 为 0-360° 服 从 均匀 分 布 。 
探测 器 参考 LHAASO-WCDA1， 各 项 参数 均 保 持 一 改 ， 在 坐标 系 中 探测 器 阵列 区 域 为 (-75 


m<x<75 m, -75 m<y<75 m), 在 数据 中 加 入 60 kHz 的 噪声 并 利用 遍历 法 取 200 ns 的 时 间 窗 口 
去 噪 之 后 信 品 比 约 为 70/22。 分 别 用 重心 法 和 椭圆 拟 合法 来 重建 事例 并 比较 其 重建 精度 和 事 
例 利用 率 。 首先 需要 验证 本 算法 对 于 蕊 位 位 于 探测 器 阵列 边沿 事例 的 重建 能 力 , 选取 真实 芯 


PE XoY) 位 于 (75 m<x<85 m,-75 m<y<75 m) 探测 器 阵列 外 紧 贴 边界 的 能 量 为 10-100TeV 
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子 事例 重建 蕊 位 , 得 到 结果 如 图 6 所 示 。 从 图 中 可 以 看 出 ,重建 蕊 位 与 真实 蕊 位 分 布 基 本 
一 致 ， 从 右 图 可 以 看 出 大 多 数 重 建 芯 位 的 偏差 dr= 4y (Xet = Xer)? F (Vet = Yir)? 小 于 10m, 
取 dr 的 68% 作 为 分 辩 率 约 为 9.6 Ko 


Core Position Ture Core Position Reconstructed 250 
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图 6: 真实 芯 位 和 重建 芯 位 分 布 《 左 ) ， 重 建 芯 位 偏差 统计 CA) 


Figure 6: Real core position and reconstruction core position distribution (left), reconstruction core position deviation statistics (right) 
为 了 定量 描述 该 方法 重建 芯 位 的 精度 和 对 探测 器 面积 的 利用 率 ， 分 别 取 4000 个 能 量 关 
1-10 TeV. 10-100 TeV. 100-1000 TeV, 1-10 PeV, 真实 芯 位 均匀 分 布 在 (-150 m<x<150 m,-150 
m<y<150 m) 区 域 的 质子 事例 ， 用 重心 法 和 椭圆 拟 合法 重建 蕊 位 。 由 于 真实 忌 位 均匀 分 布 在 
300 mx300 m 的 区 域 ， 而 探测 器 阵列 区 域 为 150 mx150 m， 因 此 取 重 建成 功 的 事例 数 与 芯 位 


位 于 探测 器 阵列 的 事例 数 之 比 作 为 探测 器 面积 利用 率 7 = a 取 dr AHI 68 多 作为 重建 芯 


位 分 辨 率 ， 得 到 结果 如 下 图 7 所 示 。 

从 图 7 中 可 以 看 出 对 于 所 有 能 量 段 的 质子 事例 ,椭圆 拟 合法 蕊 位 重建 精度 都 比重 心 法 要 
好 , 并 且 顶 圆 拟 合法 能 够 显著 提高 探测 器 阵列 利用 率 ， 并 且 随 着 能 量 增 大 其 增幅 更 加 显著 。 
对 于 1-10 TeV 的 事例 椭圆 拟 合法 的 探测 器 面积 利用 率 略 小 于 1, 这 意味 着 只 有 少数 区位 落 在 
探测 器 边沿 的 事例 重建 不 好 ， 面 积 利用 率 随 着 能 量 增 大 而 增 大 ， 对 于 能 量 为 1-10PeV 的 事 
例 其 面积 利用 率 大 于 1.3， 这 意味 着 约 有 30% 位 于 探测 器 阵列 外 部 事例 的 芯 位 重建 符合 精度 
要 求 。 左 图 中 椭圆 拟 合法 重建 芯 位 分 辨 率 随 着 能 量 增 大 而 增 大 的 现象 不 符合 预期 , 这 是 由 于 
把 大 量 芯 外 事例 纳入 统计 量 导 致 面积 利用 率 增 大 造成 的 , 通过 归 一 化 面积 利用 率 后 的 区 位 分 


状 率 可 以 得 到 证 实 。 
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图 7: 椭圆 拟 合法 和 重心 法 芯 位 重建 分 辨 率 〈 左 ) 和 探测 器 阵列 面积 利用 率 (A) 
Figure 7: Ellipse fitting method and center of gravity method core position reconstruction resolution (left) and detector array area 


utilization (right) 

各 能 段 面 积 利 用 率 都 为 0.67 的 情况 下 椭圆 法 和 重心 法 的 重建 分 辨 率 如 图 8 Atos, EE 
建 效 率 相同 的 情况 下 , 椭圆 拟 合法 的 蕊 位 分 辨 率 显 著 优 于 重心 法 ; 并 且 两 种 方法 的 蕊 位 分 辨 
率 都 随 着 能 量 增加 而 降低 , 其 中 椭圆 拟 合法 由 于 本 身 能 够 排除 大 部 分 噪声 干扰 ,其 重建 精度 
随 能 量变 化 不 显著 。 
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图 8: 面积 利用 率 为 0.69 时 各 能 段 的 芯 位 重建 分 辩 率 


Figure 8: Reconstruction resolution of core positions for each energy segment when the area utilization ratio is 0.69 
4. 结论 和 展望 

通过 蒙特 卡 罗 模 拟 数据 检验 , 椭圆 拟 合法 相 比 重心 法 具有 明显 优势 , 椭圆 拟 合法 的 优势 
主要 体现 在 两 个 方面 ， 一 是 充分 考虑 了 EAS 的 对 称 性 特征 ， 利 用 其 对 称 性 能 够 准确 重建 探 
测 器 阵列 边沿 乃至 探测 器 区 域 之 外 的 事例 , 二 是 通过 分 组 重建 的 方式 有 效 消除 了 噪声 对 已 位 
重建 的 影响 ， 从 而 显著 性 提高 世 位 重建 的 精度 。 从 150 mx150 m 的 WCDA-I 模拟 数据 的 如 
建 结果 来 看 , 椭圆 拟 合法 重建 蕊 内 事例 芯 位 重建 精度 可 以 达到 5 米 以 下 , 各 能 段 分 辨 率 不 到 
重心 法 的 50%， 重 建 精度 相 比 于 重心 法 有 显著 提高 ， 另 一 方面 能 够 解决 重心 法 无 法 准确 重 
建 探测 器 阵列 边沿 事例 蕊 位 的 问题 ， 重 建 低能 事例 的 面积 利用 率 约 比重 心 法 高 40%， 重 建 


ah 


如 
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高 能 事例 的 面积 利用 率 则 接近 重心 法 的 2 倍 ; 相 比 于 重心 法 其 蕊 位 重建 精度 和 探测 器 阵列 面 
积 利用 率 两 方面 都 有 显著 改善 。 

顶 圆 拟 合法 的 噪声 抑制 能 力 源 自 于 EAS 次 级 粒子 触发 的 Nor 信号 等 高 线 为 一 系列 同心 
椭圆 而 噪声 信号 均匀 分 布 , 可 以 根据 这 一 原理 设计 去 噪 算法 , 在 遍历 法 去 噪 的 基础 上 利用 世 
位 信息 实现 深度 去 品 。 此 外 发 现 前 锋面 次 级 粒子 到 达 时 间 也 关于 簇 芯 对 称 分布 , 可 以 基于 这 
一 特征 通过 分 区 重建 方向 规避 锥 面 修正 参数 依赖 带 来 的 问题 , 并 且 利用 方向 信息 实现 进一步 
的 深度 去 品 。 后 续 研 究 主 要 治 着 这 一 方向 展开 。 
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A New Core Position Reconstruction Method Based on Ground 


Particle Detector Array 


Abstract: Core position reconstruction is the prerequisite and basis for physical analysis of ground cosmic ray 
detection experiments. Based on the symmetry of the extended air shower and the characteristics of the detector, 
this paper proposes a new core position reconstruction algorithm—ellipse fitting method, which is used to 
reconstruct the core position of simulated data and compare it with the center of gravity method. For the event of 
cosmic ray protons in the core with an energy greater than 1 TeV, the core position resolution is less than 5 m, 
which is significantly better than the center of gravity method; for the event with an energy greater than 1 TeV 
within a range of 20 m near the array boundary, the core position resolution is less than 10 m. It has significant 
advantages over traditional methods. Because this method can reconstruct edge events and some out-of-core events 
under the premise of ensuring accuracy, the utilization rate of the events is greatly increased. For the 1-10 PeV 


proton event, the area utilization rate of the detector array is about twice that of the center of gravity method. 


Keywords: Cosmic rays; Extended air shower; Events reconstruction; symmetry; Water Cherenkov array 
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